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STAGNATION-POINT HEAT TRANSFER TO BLUNT 

SHAPES IN HYPERSONIC FLIGHT, 

INCLUDING EFFECTS OF YAW 

By A. J. Eggers, Jr., C. Frederick Hansen, 
and Bernard E. Cunningham 


SUMMARY 


An approximate theory is developed for predicting the rate of heat 
transfer to the stagnation region of "blunt Bodies in hypersonic flight. 
Attention is focused on the case where wall temperature is small compared 
to stagnation temperature. The theoretical heat-transfer rate at the 
stagnation point of a hemispherical body is found to agree with available 
experimental data. The effect of yaw on heat transfer to a cylindrical 
stagnation region is treated at some length, and it is predicted that 
large yaw should cause sizable reductions in heat-transfer rate. 


INTRODUCTION 


It has been suggested (see refs. 1 and 2) that blunting or rounding 
the leading edges of wings and bodies might substantially alleviate aero- 
dynamic heating of these regions in hypersonic flight. There is, of 
course, the added advantage that round leading edges are structurally 
more practical than sharp leading edges, especially when the problem of 
absorbing heat is considered. Another consequence of blunting may be 
increased pressure drag. In the case of ballistic vehicles, this conse- 
quence is often an advantage (see ref. l) . In the case of glide vehicles, 
however, or more generally any vehicles required to operate for sustained 
periods in more or less level hypersonic flight, increased drag may be 
viewed as a disadvantage. 

Now, to be sure, rounding or blunting the nose of a body does not 
always increase drag. Indeed, small amounts of blunting may reduce the 
drag of a body (see, e.g., refs. 3 and 4 ). The same, however, cannot be 
said for blunting the leading edge of a wing. Even small blunting causes 
a sizable increase in drag. It is natural, then, to look for methods of 
minimizing this drag penalty, and the possibility of yawing or sweeping 
the leading edge comes to mind. Impact pressures should be, according to 
simple-sweep theory, decreased in proportion to the cosine squared of the 
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angle of sweep; hence, as is intuitively obvious, large sweep should 
substantially reduce the drag penalty due to blunting. In view of this 
possibility it is important to inquire of the effect of yaw or sweep on 
heat transfer to a blunt leading edge . 

The purpose of this paper is to investigate theoretically the heat 
transfer to the stagnation regions of bodies in hypersonic flight, includ- 
ing the effects of yaw, by a simplified method which is suited to take 
account of real gas effects such as dissociation. This method, which was 
previously given limited distribution, is used along with recent estimates 
of transport properties for high temperature air, and the solutions are 
compared with some heat transfer results for blunt shapes. 


SYMBOLS 



integration constants 


specific heat at constant pressure, ft -lb /slug °R 


h 


specific enthalpy, ft-lb/slug 


k 


coefficient of thermal conductivity, ft-lb/ft-sec °R 


M 


Mach number, dimensionless 


n 


exponent of temperature in thermal conductivity and viscosity 
functions (see eqs . ( 37 ) and ( 38 )), dimensionless 


Nu 


Nusselt number based on a length 2R b and stagnation tempera- 
ture conditions, dimensionless 


P 


static pressure, lb/ft 2 (unless otherwise specified) 


Pr 


Prandtl number, dimensionless 


q heat flux per unit area, ft-lb/f t 2 -sec 

q( 0 ) heat flux per unit area at zero yaw, ft-lb/f t 2 -sec 

q(A) heat flux per unit area at yaw angle A, ft-lb/f t 2 -sec 


R 


gas constant, ft-lb/slug °R 
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R b radius of curvature of body at the stagnation point, ft 

R s radius of curvature of the shock wave at the stagnation stream- 

line, ft 

Re Reynolds number, based on twice the radius of curvature of the 

body at the stagnation point, dimensionless 

r,0,cp spherical coordinates, feet, degrees, and degrees, respectively 

T static temperature, °R 

T b temperature of the body, °R 

T 0 temperature at the interface, x = 0, with body at zero yaw, °R 

T Q (A) temperature at the interface, x = 0, with body at angle of yaw 

A, °R 

T r recovery temperature, °R 

T-j. stagnation temperature, °R 

Uqq stream velocity, ft /sec 

u,v,w velocity components in the x, y, and z directions, respectively, 
ft/sec 

u,v velocity components in the x and r directions, respectively, 
ft/sec 

x,y,z Cartesian coordinates, ft 

x,r cylindrical coordinates, ft 

5 flow deflection angle, deg 

f - R b x b - x 

€ dimensionless coordinate, — =: or — - 

R b R b 

ratio of specific heat at constant pressure to specific heat at 
constant volume, dimensionless 
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a function of density change across a shock wave, 
dimensionless 


(Ps/Poo) + 1 
(P S /PJ ~ f 


angle of yaw, deg 
density, slugs /cu ft 



a function of the coefficient of thermal conductivity 


and of temperature, ft-lb/ft-sec (unless otherwise specified) 

acute angle of shock wave relative to stream velocity vector, 
deg 

coefficient of viscosity, slugs/ft sec 

coefficient of viscosity at temperature T Q , slugs/ft sec 
coefficient of viscosity at temperature T 0 (A), slugs/ft sec 


Subscripts 

conditions just behind shock wave on the stagnation streamline 
conditions at the stagnation point of the body 

conditions at the interface between regions 1 and 2 on the stag- 
nation streamline (see sketch (b)) 


conditions in the free stream 
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Superscripts 


first derivative with respect to the x coordinate 
second derivative with respect to the x coordinate 


THEORY 

General Equations in Cartesian Coordinates 


The analysis proceeds from the equations of momentum, continuity, 
energy, and state for continuum fluid flow. The x, y, and z momentum 
equations are, respectively. 
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The continuity equation is 


I 2 + A( P u) + ^-(pv) + |-(pw) = 0 
at ax oy oz 


oo 


and the energy equation is 


Bh 


bh 


p[U b^ + V “ + W 


dy 





Sp 

Bx 
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v 


B P 
by 
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3 \Bx by Bz J 


(5) 


while the equation of state is taken in the form 


p = p(p,T) 


(6) 
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Derivations of the momentum and energy equations are given in numerous 
sources (see, e.g., refs. 6, and 7)* Note that the coefficients of 
viscosity and thermal conductivity, and the heat capacity have been 
treated as variables. It is intended that by so doing a more accurate 
solution will be obtained for hypersonic flows with their characteristi- 
cally large temperature and pressure gradients. 

Let us now consider the particular flows of interest in this paper, 
namely, those in the region of a stagnation point. 


Model of Flow and Method of Analysis 


It is instructive in setting up the model to consider the qualitative 
aspects of temperature and velocity variations in the flow along the stag- 
nation streamline. Restricting the analysis to steady hypersonic flow, 
that is H^sin 5 >>1, we will assume that the surface temperature is low 
compared to the stagnation temperature of the air. This assumption seems 
quite reasonable since practical surface materials will probably be 
destroyed if surface temperatures are allowed to approach stagnation 
temperature. It will be assumed further that the Reynolds number of the 
flow is large enough so that heat conduction and viscous shearing in the 
shock process is distinct and separate from the corresponding phenomena 
occurring in the boundary layer adjacent to the surface of the body. 
Accordingly, temperature and velocity should vary along the stagnation 
streamline similar to the manner shown in sketch (a). 



Sketch (a) 
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There is an abrupt and large increase in temperature and decrease in 
velocity of the air as it passes through the bow shock. Proceeding from 
the shock in the direction of the body, temperature continues to increase 
slowly while the velocity decreases slowly towards zero. Near the surface 
of the body, the air temperature ceases to increase and, in fact, begins 
to fall off steeply in the direction of the body temperature. The veloc- 
ity of the flow must, of course, be close to zero in this region. 

On the basis of these observations the following simplified model is 
proposed and employed throughout this study of heat transfer in a stagna- 
tion region. 



Region 2- Low- velocity, compreeeible, viscous flow 


Sketch (b) 

Since Moo is large compared to 1, Mg is substantially less than 1 and 
the detached shock wave is located a relatively short distance ahead of 
the body surface (i.e., (x s + oc-j-, ) /Ri^ < < l) • The flow between the shock 
wave and the body surface is divided into two regions. Region 1 is taken 
as a domain of essentially nonviscous, non-heat-conducting, incompressible 
flow while region 2 is taken as a domain of very low speed, but compres- 
sible, viscous, and heat-conducting flow. It is anticipated further that 
in region 2 the u and v components of velocity will be very small. The 
component of velocity w due to yaw may , of course , take on rather large 
values . 

Now it may be demonstrated with equations (l) and (2) that d p/8y 
becomes relatively independent of x along the stagnation streamline in 
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the limit as the disturbed flow extends only a short distance away from 
the body. Inasmuch as this is the type of flow of interest here, it will 
be assumed throughout this analysis that S^/dy 2 is essentially con- 
stant along the stagnation streamline between the shock and the body. 

With these assumptions, the derivative with respect to y of the 
y momentum equation yields a differential equation that becomes tractable, 
both in regions 1 and 2, when terms that vanish in the neighborhood of the 
stagnation streamline are dropped. Approximate solutions to these simpli- 
fied y momentum equations are found for the u velocity along the stag- 
nation streamline in region 1, and for the derivative of this velocity 
along the stagnation streamline in region 2 . The constants appearing in 
these solutions are determined by matching the boundary conditions at the 
shock wave and at the surface of the body, and by matching flow conditions 
at the interface. This procedure fixes the locations of the shock wave 
and interface relative to the body. 

The energy equation is simplified in an analogous manner, and solu- 
tions valid in the neighborhood of the stagnation strea m line are found 
for regions 1 and 2. The rate of heat transfer per unit area to the 
stagnation region of the body follows from the solution to the energy 
equation for region 2. 

Let us see how these thoughts apply in the case of a two-dimensional 
stagnation region. 


Heat Transfer to a Cylindrical Stagnation Region 


Zero yaw.- This problem has been treated for incompressible flow by 
Howarth (ref. 7 ) and more recently for the compressible flow by Cohen 
and Reshotko (ref. 8). One reason for re-investigating the matter here 
is to obtain compressible flow solutions which can be extended with rela- 
tive ease to the case of a yawed cylinder. In addition it was desired 
to obtain solutions which may be better suited to account for real gas 
effects, such as dissociation. 

To proceed, then, the stagnation streamlines are taken to lie in the 
x-z plane. The origin of the coordinate system is at the interface 
between regions 1 and 2, and the shock-wave and body-surface locations in 
this plane are -x g and x-^, respectively (see sketch (b)). For the case 
of zero yaw, the z component of velocity and all derivatives with 
respect to z are, of course, identically zero. 

First a solution will be found to the steady-state y momentum 
equation near the stagnation streamline in region 1. Since the flow is 
assumed incompressible and nonviscous in this region, equation (2) 
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simplifies to 


u <k+v^=-i^ (7) 

dx by P oy 

Differentiating equation (7) with respect to y there is obtained 


u + &<*! + f\ $£) + v dfv 1 dj> 

by \bx J by bx \byj by 2 P by 2 


(8) 


On the stagnation streamline v is identically zero and, therefore, bv/bx 
is also zero. In addition, the continuity equation (eq. (4)) becomes, for 
incompressible, two-dimensional flow 


8u 8v 
bx by 


(9) 


Using this information with equation (8) , one obtains 


-u 


afu fbuf 

bx 2 W 


1 8fp 
p 8y 2 


(io) 


Treating 8 2 p/8y 2 as a function of y only, and noting that equation (10) 
becomes a total differential equation along a line y = constant, yields a 
general solution for velocity along the stagnation streamline 


Cx _ -Cx 
u = Ae + Be 


where the constants A, B, and C are related by 


4ABC 2 = | — 

8y 2 


(ID 


( 12 ) 


Note that the constants may be real or imaginary, depending on the bound- 
ary conditions . 
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Now it is anticipated that the velocity u will very nearly vanish 
at the interface x = 0 (i.e., in the sense that u Q /u s < < l) ; hence B 
will be approximately -A, and the corresponding approximate solution for 
velocity is 1 


u = 2A sinh Cx 


( 13 ) 


To the same order of approximation, the second derivative of velocity at 
the interface, u 0 !1 , also vanishes. The product 2AC is just the velocity 
derivative at the interface and can be evaluated from equations (10) 
and ( 13 ), thus 


2AC = u 0 T = ± 



( 1*0 


Note that the negative root correctly describes the flow in the coordi- 
nate system of sketch (b), since velocity decreases with increasing x. 

Consider next the steady-state y momentum equation near the stagna- 
tion streamline in region 2. In this domain viscous terms must, of course, 
be retained and thus the derivative of equation (2) with respect to y 
yields 


b 2 v , bp bv t bv Su 

pu + u — — + p + 

dxdy by bx bx by 



v 

by by 




2 b 2 / du + dv \ 

3 by 2 _ ^ W + byj _ 




( 15 ) 


Now close to the surface of the body the left-hand side of this expres- 
sion is negligible and the right-hand side simplifies so that the equa- 
tion may be written (see Appendix A) 


A ( Aa = _ 

bx V bx 2 ) by 2 


(16) 


1 In the limit of zero boundary-layer thickness, this solution is 

exactly the one to which equation (ll) reduces. 
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Along the stagnation streamline this equation integrates to 




+ 


D 


(IT) 


The constant D is zero since d 2 u/dx 2 = u 0 " = 0 at the interface (x = 0) . 
Near the surface of the body, equation (17) can be integrated to obtain 




do u o 1 


(18) 


In order to satisfy the boundary condition at the body surface 
^ 1 ) = = 0 , it follows from equations (l 8 ) and (l4) that 



2q 0 



(19) 


Now p and c^p/Oy 2 can be evaluated at the shock wave since both are 
considered constant throughout region 1. In Appendix B it is demonstrated 
that for two-dimensional flow 


p 8fp \ _ _ 6 P°o 2 U co 2 

P (7 S - 1)R S 2 


( 20 ) 


where R s is the radius of curvature of the shock wave in the stagnation 
region. Substituting equation (20) in equation (19) we obtain 


fb = As-^ Vpq m i/2 2 

Rb \ 6 ) \doo Rb / R eoo 1/2 


( 21 ) 


where Re^, is the free-stream Reynolds number based on 2Rb, twice the 
radius of curvature of the body at the stagnation point. Note also that 
the effective value of 7 , the ratio of specific heats, at the shock wave 
is allowed to vary from the free-stream value. In this way, changes in 
internal molecular energy which are manifest at the high temperatures 
encountered in hypersonic flight can be considered. 
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There remains the problem of solving the energy equation. In 
region 1, the energy equation is simplified by neglecting all the viscous 
and heat-conduction terms . Then, for the two-dimensional problem con- 
sidered here, equation (5) reduces to 


u 




(22) 


for which the solution is 



(23) 


It can he seen from equation (23) that the interface temperature T 0 is 
approximately the stagnation temperature T-j- , since the velocity at the 
interface nearly vanishes. The stagnation temperature is, of course, 
given by the integral equation 



7ooRWx> 2 

2 


( 24 ) 


where for very high velocity flow the lower limit of the integral will be 
neglected. 


Next consider the energy equation in region 2. Proceeding in a man- 
ner analogous to that used in studying the y momentum equation in this 
region, we neglect the terms with the factors u, v, du/dx, and dv/dy. 
Thus equation (5) becomes simply the heat-conduction equation 



(25) 


The coefficient of thermal conductivity, k, is considered a known function 
of temperature (pressure is essentially constant) . Thus a new function of 
temperature, tj, may be defined such that 



(26) 


Then equation (25) may be expressed in terms of the function 


djj_ 

6x 2 dy 2 


0 


*n 


( 27 ) 
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Inasmuch as the body boundary is cylindrical, it is convenient to use the 
general solution to equation (27) in terms of the polar coordinates (r, 0). 
Thus 


( C n^ n + ;e) cos n0 + + ^ 

n= 1 ' 

(28) 

The origin of the coordinate system is now taken as the center of curva- 
ture of the body, and 0 as the acute angle between the radius vector r 
and the stagnation streamline. If a surface temperature is assumed inde- 
pendent of the angle 0, the solution on the stagnation streamline (0= 0) 
reduces to 


t)=A + Blnr + 


^sin n0 


n = n h + 


Blnr- + 
Rb 


C ? n 

o n r 


n=i 



(29) 


Letting — = 1 + e, where e is very small compared to unity, and 
Rb 

expanding equation (29) in a series of ascending powers of e , we obtain 


n = % 




+ 0 (e 3 ) 


(30) 


where G is the constant (B + 


2nR b n C n 


It is indicated by this 


equation that T) varies essentially linearly with e, since e 2 /2 is 
negligible compared to e and terms of higher order in e should be 
very small indeed. 3 Since e = (x b - x)/R b < < 1, equation (30) can be 
written 


n = - (b 0 ~ V 


X 

x b 


(31) 


2 The dependence of surface temperature on 9 should be small in the 
stagnation region. 

3 It should be pointed out that this argument hinges implicitly on 
the assumption that T] is a weak function of 9 near the stagnation 
streamline . 
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According to this expression, the rate of heat transfer per unit area to 
the stagnation region of the body is 


! = -V 


^o ~ % 
x b 



(32) 


The stagnation-line coordinate x-^ is substituted from equation (21) , 
and the rate of heat transfer becomes 


3 

17 Vh* 

V 


-8(7s - 1) - 

V^o 


) R b 



(33) 


A Nusselt number is defined for interface temperature conditions 
using a characteristic length equal to twice the radius of curvature of 
the body and a temperature potential of (T 0 - T^) ; thus 


Nu 


2 V R b 
* 0 (T 0 " T b) 


(34) 


or, substituting from equation (33) into (34) 


Nu 




V ) l/2 Re oo l/2 
R s/ k 0 (T 0 - Tb) 



(35) 


For a relatively cool body in hypersonic flight, it is possible to dis- 
regard the lower limit of the integral and the value of body tempera- 
ture T-^ compared to the interface temperature T Q . 

Note that the solutions given by equations (33) and (35) can be used 
for the case where viscosity, thermal conductivity, and specific heat are 
arbitrary functions of temperature. For instance, these functions can be 
calculated to include the effects of vibrational and dissociational molec- 
ular energy if the extent to which these energy modes are excited is known 
throughout the flow. It is also useful to consider the case where the 
specific heat is treated as a constant and the viscosity and thermal 
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conductivity as proportional to the nth power of temperature. In this 

/ T T 

case from equation (2^) ( note —2. « 1, — — > > 1 

\ T + 


1 ■ ©«-• 


(36) 


Noting that 


k 0 (T 0 - Tb) 


r 

I kdT 


To , N n 
T 


T 0 " Tb 


n+i 


dT = 


n + 1 


1 - ( T b /t 0 


1 " T b /T Q 


(37) 


and that 


C\ 

y-oj 


1/2 


. ,^ n/2 - m '* 


V7oo R . 


(38) 


it is seen that the expression for Nusselt number (eq. (35)) becomes 


Nu = 


1 


n + 1 Vs 


t^TJ VI 


1 / 4 / RbY /2 (l cA a/S RV /2 1 

r...E 


T b/ T o) 


n+i 


Mj 1 l-(T b /T 0 ) 


(39) 


and the rate of heat transfer per unit area to the stagnation region of 
the body is, in terms of free-stream conditions, 



3 

i/4 a, x 1 / 2 / „\S+i _ 1 / 2 ., n+2 

(S b\ Aoo R V Re oo Moo 

V Rb ) 


V R s/ V.2Cp/ n + 1 



(«*0) 


These considerations complete the zero-yaw analysis. However, before 
undertaking the study of effects of yaw on heat transfer it is appropriate 
to make a few remarks. There is the general question of the legitimacy of 
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the several assumptions underlying the present treatment of stagnation- 
point flows . In order to shed some light on this matter it is undertaken 
later in the report to examine the solutions obtained to see whether they 
are consistent with these assumptions and with pertinent results obtained 
by others . In this regard it is shown that the presumption of a constant 
second derivative of pressure normal to the stagnation streamline yields 
solutions for the distance between shock wave and body which are quite 
close to observed values. Next, it is demonstrated that, as assumed, the 
velocity u is negligibly small throughout region 2 under continuum flow 
conditions. Then it is shown that the largest of the viscous dissipation 
terms neglected in the energy equation for region 2 is indeed small com- 
pared to the heat -conduct ion terms . It is found too that the analysis 
predicts an amount of heat convected into region 2 which is the proper 
order of magnitude to account for the heat transferred to the body. 
Finally, it is shown that under comparable conditions equation (35) of 
this paper predicts essentially the same heat transfer as references 7 
and 8. 

In view of these results it would seem that the simplified analysis 
presented here for stagnation-region flows is, while on the one hand cer- 
tainly approximate, on the other hand quite capable of predicting useful 
information. Accordingly, we proceed to the study of effects of yaw on 
heat transfer. 

Yaw.- In this case the x direction is normal to and the z direc- 
tion is parallel to the stagnation line of the body (see plan view, 
sketch (c)). Then the z component of velocity has a finite value, but 
all z derivatives are again zero . 



Region I — Incompressible, nonviscous flow 
Region 2 — Low velocity , compressible, viscous flow 


Sketch (c) 
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The y momentum equation in region 1, differentiated with respect 
to y, takes the same form as equation (10) on the stagnation streamline. 
Thus the velocity u is again given by the solution 


u = 



(41) 


The z momentum equation for the stagnation streamline in region 1 
becomes, on dropping the negligible terms from equation (3 ) , 


u = 0 (42) 

dx 


which has the solution 


w 2 = 7 RT M 2 sin 2 ?\ 

oo oo oo 


(43) 


since the transverse component of velocity is unchanged on passing through 
the shock wave . 

The energy equation for the stagnation streamline in region 1 reduces 
to a form similar to equation (22) 


8u 


u ^ + w r 1 + C 


<3w 


8x 


8x 


|X = o 

OX 


(44) 


which has the solution 



u 2 + 7 00 RT og M 00 2 sin 2 A 
2 


(45) 


where again the stagnation temperature T t is given by equation (24) . 
At the interface where the velocity u is negligible , the temperature 
T 0 (A) is given by the solution to 


L 


T o 00 


p dT = — 

U P ai p 


(46) 
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which, for a constant heat capacity, Cp, is 


To(A) 


The differentiated y momentum equation for region 2 takes on the 
same form on the stagnation streamline as equation (l6) . Hence, the 
solution is 


7qqR 

2C-n 


M D0 2 cos 2 A 


(47) 




du 

Ox 


dfp xf 
dy 2 2 


o 


(A) 


1 

p 8y 2 


and the body stagnation point coordinate is 


[ x b(A )] 2 = 


2Hq(A) 

/ ^ 2 P 

J~ P £y 2 


(48) 


(49) 


Now, however, the second derivative of pressure is a function of the 
angle of yaw (see Appendix B), 



(7 S - 1)R S 2 


(50) 


so the stagnation-point coordinate is given by 


*b(M A, - 1V /4 /n o (70 R£\- /s 2 

Rb v t. V V E b ; Re„ l/2 oos 1/ ^ 


(51) 


In region 2, the solutions to the z momentum equation and the energy 
equation are considered simultaneously. The z momentum equation simpli- 
fies (to the order of this analysis), in the region of the stagnation 
streamline , to 




( 52 ) 
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Similarly, the energy equation near the stagnation streamline in region 2 
may he written (note that dw/dy is zero by symmetry) 


<5 2 t| Aw\ 2 

dx 2 Oy 2 ^ Vdx/ 


0 


( 53 ) 


In order to facilitate the solution of equation (52), it is helpful 
to observe that the yawed boundary layer, identified with the w compo- 
nent of velocity, resembles the boundary layer on a flat plate. It might 
be anticipated then that, just as in the case of the flat plate, the 
variation of w with x is relatively insensitive to variations of [± 
with x. In this event equation (52) has the approximate form 


d 2 w + d 2 w 

dx 2 + ay 2 


(54) 


The solution is taken in polar coordinates in order to conveniently fit 
the boundary condition that w is identically zero at the body surface . 
Then following the same arguments used in deriving equations (29) and (30), 
one obtains on the stagnation streamline 


w = B In 


00 

Rv, Z_, 


C n r 


•n 


n=i 



( 55 ) 


where again e = (r/Rb) - 1 < < 1. If second order and higher terms in e 
are neglected, the z component of velocity on the stagnation streamline 
becomes, in terms of x/x-^, 


w = w Q 




( 56 ) 


whence 


<3w 

dx 


Wq 

x b 


( 57 ) 
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If this result is substituted for the last term in equation ( 53 ) , 
the energy equation becomes 


d 2 q d 2 q 
dx 2 dy 2 



0 


( 58 ) 


A solution for equation (58) which satisfies symmetry conditions on the 
stagnation streamline and also the boundary conditions that tj and p are 
constant along the surface of the body is 


t](A) 


00 


% 


+ B In — + V 

Rb L, 

n=i 


C n f n 


1 




where p is a mean value of p in the stagnation region. If equa- 
tion (59) is expanded in terms of e, q takes the following form on the 
stagnation streamline 


t}(A) 





+ 0(e 3 ) 


(60) 


The constant J is evaluated by letting T| be T ) 0 when e is e Q = 
and is given by the relation 


J = 




+ ^ W ° 2 
2e c 2 


(1 + e Q ) 


+ 


(61) 


The rate of heat transfer per unit area to the stagnation region of the 
body at angle of yaw A is, from equation (60) , 


q(X) 




dx 


b 


i_ 6n(A) = i_ f j _ dw 0 2 \ 

Rb a>e b R b A 2e 0 2 ) 


(62) 


Substituting 
terms of the 


equations (26) and (6l) into this expression and neglecting 
order e Q compared to 1, one obtains 


q(A) 


1 

*b(A) 



k dT + 


M-w 0 " 


2 


( 63 ) 
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Multiplying by 


2Bh 

* 0 (T 0 - %) 


and. substituting from equation ( 51 ) yields 


gq(A)Rb 

k 0 (T 0 -T b ) 



/ ,1 R \ x 2 1/ ' 2 3-/2, 

f ^oo % ) Re m cos A 
W*) R s/ k 0 (T 0 - T b ) 



(64) 


For a constant heat capacity it follows from equations (23), (43), 
and (47) that 


T oW 

To 


wo 2 




cos 2 A 




2CpSin 2 A 

J 


(65) 


If, in addition, the thermal conductivity is proportional to the nth 
power of temperature , then 


^Tq(A) 


4 q(Tq " T b ) Jqi 


k dT = 


221 + 2 -. 
COS A 


(n + 1)(1 - Tb/T 0 ) L \T 0 cos 2 A 


1 - 


n+i -1 


(66) 


and 


P V Q 2 

2k 0 (T 0 - T b ) 


jl \ Pr sin g A 
A - T b /T 0 


p. Pr cos 2n A sin 2 A 
PoM 1 - Tb/T 0 


( 67 ) 


MCA TN 4229 


23 


Thus equation (64) becomes 


2q(A)Rb 
k 0 (T 0 " ^b) 


i / 6 \ lA ( M l/2 (SEeY 2 R e oo l/2 cos n+l/2 A 
n + 1 Vs ' V Voo R / M oo n (l - T b /T 0 ) 



1 


Tb 


T Q cos 2 A 


n+i 


+ (n + l)Pr 


M-q(A) 


sin 2 A 


} 


( 68 ) 


The ratio of equation (68) to equation (39) is the ratio of the rate 
of heat transfer to the stagnation region of a yawed body to the rate of 
heat transfer to the stagnation region of the same body at zero yaw. This 
ratio is 


n+1/2 

q(A) cos A 


<l(0) l-(Tb/T 0 ) 


n+i 


cos 2 A 


1 - 


T 0 cos A 


n+i n 


+ (n +l)Pr — 7—- sin 2 A i- 

u q (a) J 
(69) 


An analogous expression can be obtained for the ratio of Nusselt 
numbers , thus > 


Nu(A) _ ^0± f^of T q - T b 
Nu(0) q(o) \?r J T r - Tfa 


(70) 


where from equation (63) the recovery temperature, T r , is the solution to 



(71) 


However, it should be noted that the assumptions used in the analysis 
tend to be violated when the body temperature approaches recovery condi- 
tions. Therefore it should not be expected that equation (71) will yield 
accurate values for recovery temperature . 
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There remains, of course, the problem of determining p. For the 
purposes of this report p will be taken as the arithmetic average 
between p 0 (A) and p-^, that is, p = [p Q (A) + pj-J/2 . 4 In this event 

equation ( 69 ) can be written 


q(^) 

a(o) 


n+1/2.. 
cos A 


1 - (Tb/T 0 ) 


n+x 


cos 2 A 



+ 


n+ 1 
2 


Pr 


1 + 


Tv 


TqCos^A 


n 


• 2 , 
sin A 


} 


(72) 


which in the case of a relatively cool surface (i.e., T^/TqCos^ < < 1 ) 
becomes 


q.(A) 

1(0) 


n+1/2 

cos 


A 



n+ 1 
2 


Pr sin 2 A 


(73) 


Heat Transfer to an Axially Symmetric Stagnation Region 


The methods used to calculate the rate of heat transfer to a cylin- 
drical stagnation region can also be applied to the stagnation region of 
a spherical body. This analysis is parallel to that for the cylinder at 
zero yaw and thus the x axis is taken as the stagnation streamline and 
the origin of the coordinate system is placed at the interface between 
the assumed incompressible nonviscous region 1 and the viscous, low- 
velocity, compressible region 2. For the purpose of obtaining the solu- 
tions for velocity in regions 1 and 2 on the stagnation streamline, it 


^Actually this procedure might better be considered the first step in 
an iteration method where p is recalculated on the basis of the preced- 
ing calculation of T as a function of x. This refinement is not con- 
sidered warranted here where only the gross effects of yaw for angles of 
yaw well less than 90° are of principal interest. As the angle of yaw 
approaches 90 °, the analysis as a whole tends to break down due to the 
violation of the several assumptions predicated on the flow being hyper- 
sonic normal to the axis of the cylinder. 
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is most convenient to consider the momentum and continuity equations in 
cylindrical coordinates (sketch (d)). Because of axial symmetry, all 



Region I - Incompressible, nonviscous flow 
Region 2- Low-velocity, compressible, viscous flow 

Sketch (d) 

properties are independent of the angular coordinate and, accordingly, 
the r direction momentum equation becomes 


P 



+ v 


Sv 

Sr 


Sp _ 2 _S_ 
Sr 3 Sr 




2 


j£ /Sv v\ | S 

/Su Sv\ 

r \Sr tJ Sx 

_ \Sr Sx/ _ 


(7M 


While the continuity equation is 


_S^ 

Sx 


(pu) + “ (prv) 


= 0 


(75) 


In region 1 where the viscous terms are considered identically zero, 
the r momentum equation (eq. (7*0) reduces to 
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Differentiating equation (76) with respect to r and dropping terms with 
factors v and dv/dx, which vanish on the stagnation streamline, gives 



+ u 



dxdr 


1 6j? 
P dr 2 


(77) 


Now the continuity equation (eq. (75)) expands to 

P + PI + P = 0 (78) 

dx r dr 


however, on the axis of symmetry, neglecting terms higher than second 
order in r , 


dv _ v 
dr r 


(79) 


Thus, for incompressible flow, the continuity equation on the stagnation 
streamline reduces to 


du dv 
dx + 2 dr 


(80) 


Substituting equation ( 80) into equation (77) yields 


1 ( d u _ E 8 2 u _ _ 1 6 p (8l ) 

4 \dx / 2 8x 2 P dr 5 


which, upon differentiating with respect to x, and assuming 
— i-E = constant, becomes 

P dr 2 


d^u _ 

u = 0 

;w3 


(82) 


For nonzero values of the velocity u, this differential equation has as 
a solution 



+ U Q 'X + Uq 


u 


( 83 ) 
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The value of velocity at the interface, u Q , is again considered very 
small. Thus, from equation ( 8 l), the first derivative of velocity at the 
interface is, approximately. 


u . = . /I 

° P dr 2 


(84) 


As can be seen from the solution for velocity, equation ( 83 ), the second 
derivative of velocity is constant. Therefore the second derivative of 
velocity may be evaluated from equation ( 8 l) using conditions just behind 
the shock wave , thus , 


u^ 


— + —(ff) 

2u s u sPs \8r / s 


(85) 


Substituting for the values of velocity, velocity derivative, and 
second pressure derivative behind the shock wave (see Appendix B) yields 

„ „ 8 (3 - gr s ) 

0 -v 2 _ i n 2 


and 

4.7 2(7 s - d" 

U Q - - R 

7 s + 1 R s 


Now in region 2, the viscous terms are retained in equation (74) . 
Following the procedure used in studying two-dimensional flow (see Appen- 
dix A) , the r momentum equation, differentiated with respect to r, is 
simplified to 


( 86 ) 


(87) 


8 fu\ = _ 6 fp ( 88 ) 

8x \2 8x 2 J dr 2 

which integrates to 


iill = _lj?x+A 
2 6x 2 dr 2 


( 89 ) 
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and, as in the case of the two-dimensional flow. 


8u 
2 8x 


-v.2 2 

O p X 

—4 — + Ax + B 
8r 2 2 


(90) 


The constants A and B are again determined by matching the first and 
second derivatives of velocity at the interface. Thus 


A = 


^o u o" 

2 


B = 


^oV 

2 


(9D 


At the body 8u fbx vanishes, and solving for the coordinate x^ from 
equation (90) results in 


x b = 


M-o u o" 




/ 


1 + 


48 2 p 

dr 2 



In Appendix B it is shown that 



8 pJJJ 5 
r s + 1 R s 2 


(92) 


(93) 


Thus from equations (86), (87), and (93)? it can be shown that 

Uq* = (7 - l) 2 Jg(7 - 1) Rg. Reoo 
4 8r 2 p 0 u 0 " 2 (3 - 27) 2 ^o R b 


(9^) 


which is large compared to unity for any reasonably large value of 
Reynolds number (of the order of hundreds or greater). Therefore, 
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quantities of the order of unity are neglected in equation ( 92 ), the 
stagnation point coordinate reduces to 



(95) 


which is identical in form to the relation for body surface coordinate 
in the two-dimensional flow (eq. ( 19 ))* 

Next, in region 1 the viscous dissipation and heat -conduction terms 
are again neglected in the energy equation, and terms that vanish by 
reasons of symmetry along the stagnation streamline are dropped. Thus 
the energy equation for region 1 takes the same form as equation ( 23 ) for 
the two-dimensional problem and, since the interface velocity is small, 
the interface temperature T Q is again approximately the stagnation tem- 
perature T-fc . 

In region 2, the heat -conduction terms in the energy equation pre- 
dominate, and the equation reduces to the three-dimensional Laplace 
equation in the variable tj 


^2 n2 

O T| O Ti 

-I- + 

dx 2 dy 2 



(96) 


In order to fit the boundary conditions on a spherical surface, the 
solution is given in terms of spherical coordinates (?, 6, and cp) . The 
general solution which preserves symmetry about the x axis (i.e., which 
is independent of cp) is 

i = a + S + £ (c n f n + jit) f„(=°s e) (91) 

n=i 

where P n (cos 0) is the nth order Legendre polynomial in cos 0. If it 
is required that q he a constant, q b , on the surface of the body, equa- 
tion (97) can be reduced, on the stagnation streamline, to 



( 98 ) 
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then expanding in terms of e = — - 1 < < 1, results in 

-LVU 


where L is 


tj = + L(e - e 2 ) + • • • (99) 

00 

JL + \ C n (2n + l)R^ n . Neglecting the quadratic term in 

Rv, Z_i 


n=i 


€ , evaluating tj at the interface, and transforming to the variable x, 
one obtains for equation (99) on the stagnation streamline 


1 - 1 0 - ( \> - V 


( 100 ) 


Then the rate of heat transfer to the stagnation point is 


<JL=-V = 


% 


1 _ 

x b 




k dT 


( 101 ) 


which is identical in form with the zero-yaw solution for the two- 
dimensional-flow problem (eq. (3l)). Note that in Appendix B the second 
derivative of pressure given by equation (Bl8) is larger by a factor 
of 4/3 than it is for the corresponding two-dimensional-flow case with 
the same shock-wave curvature (eq. (B17))- Thus xq> given by equa- 
tion (95) is changed by the factor (3/4 ) 1 4 and the rate of heat trans- 
fer to an axially symmetric region becomes 




2 (7 B - DJ 


l/4 


Ho R s . 


1/2 Re 1/2 
00 

R-K 


k dT 


( 102 ) 


while the corresponding expression for Nusselt number is 


Nu = 


8 


la - I 


1/4 


Re, 


1/2 


Ho v k 0 (T 0 - T b ) 


k dT 


T> 


(103) 
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Examination of Analysis and Assumptions^ 


A number of assumptions have been made in the theoretical analysis, 
and it is desirable now to show that the solutions obtained are both 
realistic and consistent with these assumptions. In particular, it will 
be shown that the presumption of a constant second derivative of pressure 
normal to the stagnation streamline yields solutions for the distance 
between shock wave and body which are reasonably close to observed values. 
Secondly, it will be demonstrated that the u velocity throughout 
region 2 is indeed small, as assumed in the analysis, if the Reynolds num- 
ber is large enough for continuum flow conditions. In addition, it will 
be shown that for region 2 the viscous-dissipation terms due to the 
u and v component velocity derivatives are small compared to the heat- 
conduction terms in the energy equation, again provided the Reynolds num- 
ber is not too small. These findings, then, help to justify the manner 
in which the momentum and energy equations were treated in the analysis . 

Now it is obvious that the assumption of an abrupt transition from 
nonviscous, convective flow to viscous, conductive flow is a substantial 
idealization of the actual flow. 5 It is possible, however, to make a 
gross check on the self-consistency of this model by comparing the amount 
of heat convected across the interface with the amount conducted to the 
body surface. When this is done it is found that from a heat-flow point 
of view, the model is self-consistent (i.e., heat convected provides for 
heat conducted) . 

As a final point, a comparison will be made between the analysis of 
this paper and the heat-transfer solutions for low-velocity flow given by 
Howarth (ref. 7) and Cohen and Reshotko (ref. 8). 

Distance between shock wave and body .- Consider first axially 
symmetric flow. The velocity in region 1 was found to be (eq. (83)) 

U " 

u = u Q + u 0 T x + -75- x 2 (104) 


5 Strictly speaking, this idealized model should be considered simply 
a first approximation to the correct situation. A second approximation 
would be to divide the domain between the body and shock wave into three 
regions rather than two as was done here. 
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Then, the shock-wave coordinate must be 


u ' - u ' 


Xc = 


U r 


(105) 


It can be shown from the relations in Appendix B and equations (84) 
and (85) that 


u s ' = 


4U„ 


(7 S + 1)R £ 


u r 


U r 


. /_ it 
J P dr 2 

8(3 - 2 7) U, 


4u 


(/s + 1)R £ 


J ' 2 (7 S - 1) > (106) 


7 2 - 1 R; 


CO 

2 


Substituting these relations into equation (105) yields 


R £ 


(7 s - l)[l - J2(7 S - 1 )] 
2(3 - 27 s ) 


(107) 


Note that for 7 S = 1.5, Uq" vanishes and the velocity profile becomes 
linear. For this case x s /R s reduces to (/ s - l)/4. 

The actual distance between the body and the shock wave is, of course, 
the sum of x s and X5. However, it can be shown from equations (95) and 
(107) that Xb is small compared to x s for reasonably large Reynolds 
numbers, and x-^ will therefore be neglected. The ratio x s /R s calcu- 
lated from equation (107) for 7 S equal 1.4 is 0.105* Measurements of 
x s /Rb taken from spark photographs of high-velocity spheres presented 
by Charters and Thomas (ref. 9) and Dugundji (ref. 10) approach this 
value closely at high Mach numbers (i.e., x s /Rb about 0.11 at Mach num- 
ber 4). Heybey (ref. 11) has developed a theory which fits the data of 
references 9 and 10 closely and, for the limit of infinite Mach number, 
predicts x s /R b about 0.12. Thus it is seen that at high Mach numbers, 

the assumption that the second derivative of pressure is constant and 
that the ratio Rb/Rs is near unity yields results which are consistent 
with experimentally observed distances between the shock wave and a 
spherical body, as well as with the theory of Heybey . 
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It is of interest to calculate the shock-wave coordinate for two- 
dimensional flow as well. Recall that the solution for velocity in 
region 1 for this case is 


u 


u o 

— sinh Cx 
C 


(108) 


and thus 


— = u Q ' cosh Cx (109) 

8x 


The velocity derivatives at the shock wave and at the interface, given 
in Appendix B, are, respectively, 


u 


I 


s 


u 


. 2U 0P 

(7 a + 1 ) R s 





> 

Uco J6(y s - D 

(7g + l)^g 


Then the product Cx s is given by 

2 

Cxo = arc cosh — ~ ■ ■ 

V 6 (r s - 1) 


( 110 ) 


( 111 ) 


With Cx s known and the velocity at the shock wave 


u s 


7s - 
7 S + 1 


( 112 ) 


The shock-wave coordinate becomes 


Xc 


- 1 


Cx £ 
sinh Cx, 


s 


(113) 
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For a 7 S of 1.4, Cx s takes the value 0.75 and since the sinh func- 
tion is very nearly linear over this range, rather close hounds on the 
shock-wave coordinate are imposed by 


u s 

< 

x s 

< 


U S ’ 

V 


0 



or 


?s 


- 1 


2 




(114) 


(115) 


The exact theoretical solution for x g /R s at / g = 1.4 is 0.236. 
According to the theory then, a shock wave with given radius of curvature 
should be detached from a cylindrical body about twice as far as from a 
sphere, assuming R s /Rb 85 1. 

Magnitude of velocity in region 2 .- The y momentum equation in 
region 2 was reduced to 



b 2 p xf 
Sy 2 2 




(116) 


The left side of this equation may be approximated by — (pu) with the 

ox 

presumption that velocity in region 2 is small. Then equation (ll6) may 
be integrated to 


pu 


a % 

ay 2 6 


t^o u o' x 


Po u o 


(117) 


Solving for u Q , noting that velocity vanishes at and substituting 

from equations (l4) and (19) > one obtains 


u 


o 


- - u, 


o’ x b = 3 



2M-q u q ' 

Ps 


(118) 
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It follows that the ratio of interface velocity to the velocity at the 
shock wave is given by 



8 u q _ 1 6 / ^o\ u o * 

9 P s u s ~ ~ ~9 \iw 


(119) 


which on substituting the relations given in Appendix B for u Q ' and u s 
becomes 


fio\* = 16 / 6 /t^ RjA j_ 

V U s/ 9 7 7 S - 1 Voo R S/ Re oo 


( 120 ) 


It can be seen that for large Reynolds numbers, of the order of 
hundreds or greater, the velocity at the interface is small compared to 
the velocity at the shock wave . Since the velocity in region 2 is every- 
where less than at the interface (see eq. (ll7)), the solutions obtained 
for velocity are consistent with the assumption that velocity is small 
throughout region 2. 


Viscous dissipation in region 2 .- Although the derivative of velocity 
vanishes at the body surface, it increases parabolically (see eq. (ll6)) 
to Uq * at the interface. Since viscous dissipation terms due to this 
velocity shear were neglected in solving the energy equation, it will be 
shown that the maximum value of these terms, which occurs at the interface, 
is small compared to the heat-conduction terms like d 2 q/8x 2 (note that by 
continuity 8v/8y contributes a dissipation term of the same magnitude as 
8u/8x) . From equation (30) it can be seen that the term 8 2 tj/8x 2 is 
nearly constant everywhere along the stagnation streamline in region 2. 

Then the ratio of differential terms in the energy equation is, by 
equations (30) and (ll8) , 


4p 0 (uq') 2 

(8V8x 2 ) b 



( 121 ) 


If equation (l2l) is evaluated for constant heat capacity and thermal 
conductivity proportional to the nth power of temperature, there is 
obtained 


4m- 0 (u 0 1 ) 2 

(8 2 T)/8x\ 



= l8(n + 1) 
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Substituting for velocity ratio Uq/U^ from equations (120 ) and (b 4) 
and for the ratio R^/x^ from equation (21) , there results 


V)(u q ') 2 
(^ti/Sx 2 )^ 


= l6(n + l)Pr 



(- 


6 y /4 

7 - 1) 




-1/2 


(123) 


Once again the square root of Reynolds number is the predominant term for 
conditions of continuum flow and thus the viscous dissipation terms in 
the energy equation are small compared to the conduction terms in region 2. 

Heat convection across the interface, x = 0 .- Next consider the ratio 

of the heat convected across the interface, pu 0 j° CpdT, to the heat 

'-'o 

transfer at the stagnation point of the body, The value of u Q 

given by equation (ll8) and -% * from equation (32) yields 


PUq 

V 



(124) 


Again evaluating for constant heat capacity and the nth power tempera- 
ture function for thermal conductivity, and noting from equations (l4) 
and (19) that u^x^ 2 reduces to -(2p 0 /p), one obtains 


pUpCpTo 

V 



n + 


1) 


(125) 


This ratio is the order of unity , and thus the right magnitude of heat 
is convected across the interface to balance the heat conducted to the 
body. The above result also provides a check on the value of x^ which 
was obtained by matching u Q ' as a boundary condition of the y momen- 
tum equation. 
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Low- velocity heat transfer .- For hypersonic velocities it was found 
that taking shock-wave curvature equal to body curvature on the stagnation 
streamline gave approximately the correct answer for the distance between 
the body and the shock wave, so presumably the ratio R b /R s should be 

taken near unity when calculating the heat transfer as well. Undoubtedly 
this ratio will be somewhat less than unity for low Mach number supersonic 
flow, and it is of interest to see what the solutions developed in this 
paper will predict for this case (even though the assumptions made in the 
analysis are not expected to hold as well for the low- velocity flow con- 
ditions) . For this purpose it is convenient to express the body coordinate 
x b i- n terms of (8v/8y) 0 which by continuity equals -u Q ' . From equa- 
tions (9), (14), and (19) 


x b 


2p c 


p(8v/8y) 


(126) 


then solving for Nusselt number from equations (32) and (34) for the case 
of the cool wall (T b /T Q < < 1, and n = 1/2) one obtains 


Nu .|^.°.47B bi 


/p(dv/dy) 


(127) 


The method of boundary-layer solution for low-velocity flow about a cylin- 
der given in reference "J , yields for the derivative of velocity component 
normal to the stagnation streamline at the edge of the boundary layer, in 
the notation of this paper , 



( 128 ) 


Substituting in equation (127) results in 


Nu = 0.92 Re s 1/2 

where the small differences between p s and p Q are neglected. The con- 
stant 0.92 compares favorably with the value 0.95 given by Howarth for 
Pr = 0.72. This agreement is especially remarkable in the light of the 
fact that the analysis of reference 7 is for constant thermal properties, 
while variation in thermal properties is an essential feature of this 
analysis . 
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Cohen and Reshotko (ref. 8) find that the solution for a compressible 
boundary layer gives the following relation at the stagnation point of an 
axially symmetric body 


. 0.440 IHZ 
Db V ^ 


(129) 


for the case of a cool wall and a Prandtl number 0-7- If the radial com- 
ponent of velocity v is taken proportional to y, the ordinate can be 
eliminated and equation (129) reduces to 


Nu = 0.440 Db 


p(8v/8y) 


( 130 ) 


The factor 0.440 given by Cohen and Reshotko compares favorably with the 
factor 0.47 given in equation (127) • 


HEAT-TRANSFER RESULTS FOR BLUNT SHAPES IN HYPERSONIC FLIGHT 


Temperatures in the disturbed flow about vehicles in hypersonic 
flight may be sufficiently large to dissociate air molecules into atoms 
or even to ionize the atoms* At present the chemical reaction rates for 
these processes are not known with certainty. Available experimental 
evidence (ref. 12) indicates that air will be in equilibrium throughout 
the stagnation region flow for vehicles in flight at velocities up to 
26,000 feet per second, and at altitudes up to about 200,000 feet. At 
much greater altitudes, the atmosphere is so rarefied that the chemical 
reactions will probably be frozen and the air will behave essentially as 
a gas with constant specific heat. These two limiting cases, at least, 
can be treated within the framework of the present analytical results . 

For this purpose it will be convenient to consider the heat-transfer rate 
expressed in the form of a parameter which is relatively independent of 
scale size or density. From equation (102), such a parameter is given 
by 


. (!• . p° k dT 

sJReoo/2 VPco J VjW Jo 


(13D 


where it has been assumed that the surface temperature is negligible 
compared to T 0 and that the shock-wave curvature equals the body 
curvature in the stagnation region. Equation (131) applies in the 


NACA TN 4229 


39 


spherical case; the rate of heat transfer to a cylindrical stagnation 
region may be smaller by the factor (3/4) 1/4 according to equation (33). 
Note that the integral may be evaluated with the thermal conductivity 
coefficients taken at constant pressure, since the pressure is relatively 
invariant along a stagnation streamline . The integrals have been calcu- 
lated graphically using the data given in reference 13 and the results 
are shown in figure 1. 

For the case where chemical reactions are frozen, all translational, 
rotational, and vibrational modes of energy are considered fully excited, 
C p /r is taken a constant at 9/2, and the coefficients of viscosity and 
thermal conductivity are taken proportional to the half power of tempera- 
ture. The heat-transfer parameter given by equation (l3l) for these con- 
ditions is shown in figure 2 for flight velocities from 5,000 to 30,000 
feet per second. 

For the case of chemical equilibrium, Feldman (ref. l4) has calculated 
the densities and stagnation temperatures which occur behind shock waves, 
and reference 13 gives values for the coefficients of viscosity and thermal 
conductivity. The chemical reactions, which keep the flow in equilibrium, 
cause the thermal conductivity to be much larger than in the frozen flow, 
but this effect is compensated for by the large decrease in stagnation 
temperature due to the strong heat sinks created by the reactions. Inci- 
dentally Kuo (ref. 15) finds similar compensation for the case of heat 
transmitted through the boundary layer along a flat plate. Because of 
the compensating effects, it is not immediately apparent whether the 
integral in equation ( 131 ) will be increased or decreased by the dissocia- 
tion and ionization reactions . In all the cases calculated it is found 
that the integral is slightly greater under equilibrium conditions . In 
addition, both of the other factors in equation ( 131 ) are increased slightly 
by the chemical reactions leading to equilibrium. The density ratio across 
a normal shock may increase more than a factor of 2 (see ref. l4) , but the 
heat-transfer rate varies only as the fourth root of this ratio and is not 
strongly influenced. Reference 13 finds that the coefficient of viscosity 
is increased somewhat at equilibrium, but this also is compensated by the 
decrease in stagnation temperature. The resulting ratio |-Loo/po is 
increased slightly, but again the effect on heat transfer is weakened by 
the square-root dependence on this factor. The total result of increases 
in all factors is that the parameter qR-^Re qq’ 1 ^ 2 is the order of 30 per- 
cent greater for stagnation region flow in equilibrium than for such flow 
in which the chemical reactions are frozen. The difference is indicated 
by the two curves in figure 2. 

The heat transfer calculated for the equilibrium flow is in satis- 
factory agreement with the experimental results reported by Rose and 
Riddell (ref. l6) as indicated in figure 2. It may be noted that there 
is a few percent change in the heat-transfer parameter due to different 
ambient temperature and pressure conditions which occur at different 
altitudes , but in view of the order of the approximations inherent in the 
theory and of the ±20-percent variation in experimental results, the 
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change is not significant enough to be shown in figure 2. The theoretical 
results for. equilibrium flow also agree with numerical integrations of 
more complete boundary- layer equations, including chemical reaction terms, 
which have been made by Fay and Riddell (ref. 17). Thus it is concluded 
that the approximate theory presented in this report retains the essential 
relationships which influence stagnation-region heat transfer. 

In view of the foregoing results, it seems reasonable that the present 
theory would also yield approximately correct values for the effects of 
yaw. Figure 3 shows the product of the secant of the yaw angle and the 
ratio of the stagnation- region heat flux at yaw to the flux at zero yaw, 
for the case where the wall temperature is negligible compared to the 
stagnation temperature. This quantity, q(A)/q(0)cos A, equals the ratio 
of the heat flux per unit of span normal to the stream velocity, to the 
same heat flux at zero yaw. The ratio of the heat flux per unit area is 
just q.(A)/q(0), of course. The frozen flow case was calculated from 
equation (73) where the Prandtl number was taken equal to 0. 75, and this 
result is independent of velocity. The equilibrium flow heat transfer 
was calculated for flight at 26,000 feet per second at 100,000 and 150,000 
feet altitude from the relation 


< l (>0 

q(0 jcos A 


( 132 ) 

which is derived from equations (33) and (64). At small angles of yaw, 
the effect of yaw is to reduce heat flux slightly more in the chemically 
frozen flow than in the two equilibrium flow cases shown. This is due 
primarily to particular variations in the integral of thermal conductivity 
with stagnation temperature in the equilibrium flow (fig. 1) and is not 
necessarily typical. At larger angles of yaw, the reduction in heat 
transfer is about the same in either case. As shown in figure 3, the 
stagnation- region heat flux per unit span is reduced approximately by the 
factor (cos A) 1/2 at large angles of yaw up to 70°. The corresponding 
heat flux per unit area is reduced by about the factor (cos A) 3/2 . 

The effect of wall temperature on the reduction in heat flux caused 
by yaw is shown in figure 4. The heat-transfer rates are graphed for 
wall temperature to stagnation temperature ratios of 0.2, 0.1, 0.05, 0.02, 
0.01, and 0 for flow in which the Prandtl number is equal to 0.75 and the 
dissociation reactions are frozen in a state of no dissociation (note 
that vibrational energy may be excited, however, without appreciable 
influence on the ratio, q(A)/q(0)). At high yaw angles, the viscous 
crossflow is the predominant factor contributing to the stagnation-region 
heat transfer. The principal effect of high wall temperature is to main- 
tain sizable air temperature, and therefore sizable viscosity and viscous 
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dissipation (see eq. (58)) throughout the crossflow boundary layer. 

As a consequence, the stagnation-region heat flux per unit span does not 
decrease monotonieally with increasing yaw angle, but goes through minima 
as shown in figure 4 . As the wall temperature is reduced, the viscosity 
near the body gradually becomes negligible compared to the viscosity near 
the edge of the boundary layer (i.e., at the interface x = 0) . The 
results are not strongly influenced until the wall temperature is depressed 
to the order of 0.1 the stagnation temperature. Then as wall temperature 
is further decreased, the heat flux rapidly approaches the limiting value 
given by equation ( 73 ) • Because of strong compensating effects, similar 
to those which occur in the cold-wall case at zero yaw, it is likely that 
the effect of wall temperature on heat transfer to yawed shapes in equi- 
librium flow will be quite similar to that shown in figure 4. 


CONCLUDING REMARKS 


The theory for heat flux to the stagnation region of blunt axially 
symmetric shapes in hypersonic flight, which is developed in this report, 
is found to agree favorably with other theoretical results and with avail- 
able experimental evidence. It is concluded that this theory, though 
approximate, preserves the essential functional relationships which influ- 
ence stagnation- region heat transfer. A similar analysis is made for the 
heat flux to a cylindrical stagnation region at angle of yaw. It is 
deduced that wing sweepback should reduce the heat flux per unit area at 
the leading edge approximately by the factor (cos A) 3/2 , if the wall tem- 
perature is held relatively cool. This will reduce the cooling required 
to alleviate hot spots and the thermal- stress concentrations induced by 
heating in the stagnation region at very high-speed flight. The total 
stagnation-region cooling required for a given wing span will also be 
reduced in this case, since the heat flux per unit span decreases approxi- 
mately as (cos A) 1 / 2 . 

Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., May 2, 1955 
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APPENDIX A 


SIMPLIFICATION OF THE y MOMENTUM EQUATION IN REGION 2 


The steady- state , two-dimensional y momentum equation (see eq. (2)), 
differentiated with respect to y, yields 


d 2 v dp dv dv du 

pu + u — - — + p + 

dxdy dy dx dx dy 



dp dv 
dy dy 







(Al) 


Now on the stagnation streamline the velocity v is identically zero 
and therefore all x derivatives of v are zero. Also, all odd order 
y derivatives of functions like density p, viscosity p, pressure p, 
and velocity u vanish since, by symmetry, these functions are even. 

In addition, it is assumed that near the stagnation streamline the veloc- 
ity u is so small throughout region 2 that terms with this factor may 
be neglected. With this assumption an additional useful relation can be 
deduced from the continuity equation 


P 


du 

dx 


+ u 





(A2) 


Eliminating the terms with factors u, v, or dp /dy from equation (A2) 
there results, as for incompressible flow, 


du dv 
dx dy 


(A3) 


Note that all derivatives of the sum du/dx + dv/dy are also zero in the 
regions where equation (A3) will hold. 

Applying the above considerations simplifies equation (Al) to 
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Now it will be assumed, as is usual, that the viscous flow in the region 
of the stagnation point of a blunt body is similar to viscous flow at 
the stagnation point of a body with infinite radius of curvature insofar 
as the velocity derivatives are concerned (i.e., the principle effect of 
the body curvature is to determine the magnitude of the pressure deriva- 
tives) • Accordingly, dA/dy 3 and dA/dy 2 will be supposed to vanish 
in the stagnation region. Then expansion of the second member of the 
right side of equation (a 4) yields 

/d 2 p Sv p dp dfv 

\dy 5 dy dy dy 2 



in which the only term retained is 2(d 2 p/dy 2 ) (dv/dy) . Similar expansion 
of the last member of equation (a 4) gives 


d 2 d /du + dv\ + dp ( d 2 u + dfv \ dp /d^ + d 2 v V / d 3 u 5 3 v \ 
dxdy \dy dx/ dy Vdxdy d x 2 / dx Vdy 2 dxdy / \dy 2 dx dydx 2 / 


Note that from equation (A3), d 3 u/dy 2 dx is equivalent to -(dA/dy 3 ) and 
will therefore be neglected. The terms retained in this equation, then, 
are (dp/dx) (d 2 v/dxdy) + p(d 3 v/dydx 2 ) . These terms can be combined into 



Equation (a 4) thus is reduced to 


_d_ f djA _ _ d^p d 2 !! dv /dv\ 2 
dx V dx 2 / dy 2 dy 2 dy P \dy/ 


(A5) 


The derivative dv/dy vanishes at the surface of the body, so that in 
the immediate region of the stagnation point, equation (A5) takes on the 
approximate from 


A. L SsV - — 

dx V dx J dy 2 


(A6) 


This expression will be taken to hold near the stagnation streamline 
throughout region 2. 
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APPENDIX B 

BOUNDARY VELOCITIES AND PRESSURE DERIVATIVES 


For hypersonic Mach numbers, the density ratio across an oblique 
shock wave is 



7 S + 1 

7 - 1 

s 


(Bl) 


then the pressure just downstream of the shock is 


2p U rn 2 cos 2 a 

p s = 00 00 

7 S + 1 

where a is the acute angle between the shock wave and the normal to the 
free-stream velocity vector (see ref. l8) . It can also be shown that the 
v component of velocity just downstream of the shock is 



v 


s 


2 

7 S + 1 


btosin o cos a 


(B3) 


while the u component on the stagnation streamline is 


u 


s 


7s - 1 

7 S + 1 


U„ 


(b4) 


In evaluating the derivatives, consider a shock wave with radius of 
curvature R s . Let s be the distance along this profile measured from 
the stagnation streamline and x(s) and y(s) be the equations for the 
shock-wave coordinates. Then 


while 


dv dv dx dv dy 

ds dx ds dy ds 


(B5) 
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In terms of the radius of curvature R s , the differential equations for 
x(s) and y(s) are 





(B7) 


J 

and at the stagnation streamline (ds = 0) the following conditions hold 


dy 

ds 


1 


dx 

ds 



0 
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d 2 x _ 1 
ds 2 Rs 
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J 


(b8) 


Then, at the stagnation streamline, equations (B5) and (b 6) become 


and 



iV\ = <3. 2 p 

.^y 2 / ds 2 



Now by continuity and equations (B3) and (B9) 


(B9) 
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for two-dimensional flow. For axially symmetric flow the corresponding 
relation is 


^ = -2^ = 




\Sr, 


4U 


(7 S + l)Rc 


(B12) 


According to equation (B2) the first right-hand term of equation (BIO) is 


d^ _ _ 4p oo u oo 2 
ds2 (7 + 1)R S 2 


(B13) 


while the next term, - — ( ^ ) , is evaluated using the x momentum 

R s Vdx/ S 

equation (eq. (1)) which for the nonviscous incompressible flow region on 
the stagnation streamline reduces to 


dp du 

dx “ " PU dx 


(Bl4) 


According to equations (Bll) and (B12) , equation (Bl4) becomes 


dp ^(^s ~ 
dx (7 + 1) ; 


R c 


(B15) 


and 


dp = 4(y s - 1) p s U o 
dx (7 + l) 2 Rs 


(Bl6) 


for the two-dimensional and the axially symmetric flow cases, respectively. 
Then the corresponding second partial derivatives of pressure are 


d 2 p _ 6(7 s - l)p s U a 


dy* 


( 7 , + 1) Rs' 


(B17) 
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and 



8(7 a - DpkUJ 5 

(7 S + 1) V 


(B18) 


Note that 7 S can have values somewhat different than 1.4 if vibra- 
tional and dissociational energies are excited at the shock wave. The 
results of this appendix are consistent if y Q is defined by equation (Bl) 
from the ratio of densities across the shock wave . When additional energy 
modes are excited at the shock wave, this effective value of 7 is not 
exactly the ratio of specific heats . 

It can be seen that for the case of a yawed two-dimensional body, the 
same relations hold as for the body at zero yaw except that the velocity 
U OT is replaced by the normal component of velocity, U^cos A. Thus the 
yawed two-dimensional body has a second derivative of pressure 


d 2 p = 6 ( 7 q - l)p s U 00 2 cos 2 A 

a y 2 (r s + i) 2 Rg 2 


(B19) 


In the above relations the radius of curvature of the shock wave R s 
is yet undetermined. In the limit of infinite free-stream Mach number, 
the ratio of shock wave to body curvature, R s /Rfc>> might be expected to 
approach unity as an upper bound. On the other hand, a value of R s /Rb 
consistent with incompressible boundary-layer solutions may be a reason- 
able lower bound. In this regard Howarth (ref. 7 ) reports that for two- 
dimensional flow 


dv _ ^s 
dy R b 


(B20) 


which, according to equations (b 4) and (Bll), corresponds to a ratio 


R.= 


Rb 7s ~ 1 

Sibulkin (ref. 19 ) , using a similar analysis finds that 

bv _ 3U' S 

dy ~ 2R^ 


(B21) 


(B22) 
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for axially symmetric flow. This corresponds to the ratio 
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Figure 1.- Integral of thermal conductivity as a function of temperature, 
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Figure 2.- Heat-transfer rate to the stagnation region of a blunt, axially 

symmetric shape in hypersonic flight . 
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Figure 3>- Effect of yaw on heat-transfer rates to the stagnation region 
of a blunt cylindrical shape with a cool wall. 
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Figure 4.- Influence of wall temperature on the reduction in stagnation 
region heat flux due to yaw; chemical reactions frozen. 
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